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Abstract 

We have examined holographic renormalization group {RG) flows of the shear vis- 
cosities in anisotropic holographic superfluids via their gravity duals, Einstein-£'f7(2) 
Yang-Mills system. In anisotropic phase, below the critical temperature T c , the SO(3) 
isometry(spatial rotation) in the dual gravity system is broken down to the residual 
SO(2). The shear viscosities in the symmetry broken directions of the conformal fluids 
defined on AdS boundary present non-universal values which depend on the chemical 
potential \x and temperature T of the system and also satisfy non-trivial holographic 
i?G-flow equations. The shear viscosities flow down to the specific values in IR region, 
in fact which are given by the ratios of the metric components in the symmetry unbroken 
direction to those in the broken directions, evaluated at the black brane horizon in the 
dual gravity system. 
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1 Introduction 



The Gauge/gravity duality has shed light on strongly coupled field theories. Especially, 
Fluid/gravity duality is widely studied and it provides much useful information about con- 
formal fluid dynamics in the effectively long wavelength limit. Many crucial quantities in real 
time formalism in conformal fluids can have been obtained via retarded(advanced) Green's 
function(s) from their gravity duals [JJ [2]. 

The most celebrated example from holographic fluid dynamics is known as the ratio of the 
shear viscosity n to entropy density s. The ratio seems to be universal for many conformal 
fluids, f = 4^ [31 El El [6] , °f which holographic duals are Einstein gravities in asymptotically 
AdS space. So far, the only violation of this universality appears when string effects or quantum 
effects are taken into account in the dual gravity systemJTJ [HI E] - 

One interesting direction to investigate the shear viscosities in boundary conformal fluids 
is studying their holographic renormalization group(-RG) flow equations. It has been argued 
that Wilsonian RG flow is consistent with holographic RG in [TO] [TTj . In the dual gravity 
system, the radial direction, r of AdS space is identified with Wilsonian i?G-direction and 
that radial direction is related to the energy scale of the dual fluids. AdS boundary is treated 
as f/y-region whereas the black brane horizon is as Ji?-region. In between, one can define a 
fluid dynamics at intermediate energy scale and such conformal fluids are defined on a hyper 
surface located at r = r, where < f < oo, is the black brane horizon and AdS boundary 
is located at r = oo. 

In principle, holographic RG flow and the radial evolution of bulk Einstein equations of 
motion are different. In []JJ [TTJ [T2], however, the authors have provided a general proof that 
two flows are indeed equivalent in the limit of classical gravities. Therefore, once one solves 
the bulk equations of motion, the corresponding flows are completely known. 

For the cases so far[6l[13], the shear viscosity does not run along the radial direction. From 
explicit calculation, it is obtained that df-n = for these cases. In [14J, the shear viscosity 
runs, but the entropy density also runs in such a way that the shear viscosity to the entropy 
density ratio does not change. In [6] , the authors argue that the physical reason why the shear 
viscosity shows trivial flow may be found in membrane paradigm. The membrane paradigm 
may provide an argument that the linear response of the boundary fluid dynamics is completely 
captured by that of the horizon fluid dynamics in the small frequency limit. The transport 
coefficients on the boundary fluids can be expressed in terms of those on the horizon fluids 
only. The shear viscosity is expected to be an apparent example of this. 

However, we will provide an example of non-trivial RG flows of the shear viscosities from 
anisotropic conformal fluid dynamics which has been suggested recently [T5| IT6] . This fluid 
system displays an order parameter which depends on a certain spatial direction below its 
critical temperature T c and becomes p wave superfluid, in which the ratio in that direction 
does not show the universal value (TTJ [TB] H. The gravity dual of this conformal fluids is Einstein- 
SU(2) Yang- Mills system! It is an exact solution of the bulk action with an AdS black brane 

2 The universality of the shear viscosity and the entropy density ratio is violated for some other cases too, 
as higher derivative gravity theories [211 [25j [26l EB EEJ [29l [30] , Einstein- Axion-dilaton system [31] and a study 
on non-trivial temperature dependence of the ratio [32 

3 See [201 [23 1221 123] for pioneering works on connection between Einstein-S'[/(2)-Yang-Mills and p-wave 
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in 5-dimension with non-zero chemical potential turned on for the temporal component of 
the gauge potential proportional to 03 of SU{2) gauge group. The boundary metric enjoys 
SO (3) global symmetry (as the spatial rotation in x,y and z directions in our coordinate 
system). At high temperature (equivalently small chemical potential /j), the system stays in 
the isometric phase. However at a certain chemical potential \i = fi c (equivalently, a certain 
temperature T = T c ), 5*0(3) symmetry is spontaneously broken into residual SO (2) (rotations 
in y and z directions) because one of the spatial components^ directional component) of the 
Yang-Mills fields develops a non-trivial zero mode (as a solution of the linearized Yang-Mills 
field equations) and that mode is thermodynamically more favorable than the trivial zero 
mode(isotropic phase). It turns out that in the region of fj, > fi c , this mode condenses and 
there is a new anisotropic superfluid phase. 

In this note, we have obtained flow equations of the shear viscosities in symmetry broken 
direction in the anisotropic phase using the radial flows of the bulk equations of motion. It turns 
out that the flow equations are not trivial, since they contain contributions from interactions 
with perturbative Yang-Mills fields. 

The shear viscosities will non-trivially flow into some specific values deeply in IR region, 
which would correspond to the black brane horizon in the dual gravity system. It turns out 
that the shear viscosities at the black brane horizon are expressed in terms of the ratios of the 
metric components in the symmetry unbroken directions to the those in the broken directions. 
For our case, the metric of the bulk spacetime is still diagonal in anisotropic phase. Therefore, 
for example, the shear viscosity r} xy will be given by 

(1) 

at the horizon 

where k 5 is the 5-dimensional gravity constant, G xx and G yy are diagonal metric components 
of the bulk spacetime and is the universal value of the shear viscosity. 

The reason for this behaviors in IR region is due to the causal boundary condition at the 
black brane horizon. Near horizon, the fast oscillating factor, (r — r^) - *^ of the gravitational 
perturbations which contribute to the shear viscosity calculations has dominant contribution 
and the other regular factors are relatively suppressed, where (5 is a positive real number and 
v is frequency of the fields. In the equation of motion of the gravitational perturbations, 
the interaction terms from Yang-Mills fields are relatively suppressed near horizon and the 
equation becomes the same form of the near horizon limit of that in isotropic phase. Since /3 
depends only on the background metric evaluated at the horizon @, the shear viscosities will 
do as such. In some literatures (e.g see[3T]), the form of the shear viscosity as EqflT]) at the 
horizon is expected but they are the cases for trivial RG flows. We have shown that this is 
true for the case of the non-trivial flows too. 

In the last section in this note, we will provide a specific example for this argument. Us- 
ing the analytic solutions in Einstein-5'J7(2) Yang-Mills obtained in [19] . we explicitly show 
that Eq([T]) is valid near critical point. This analytic solution can be obtained only when the 

holographic superfluids. 

4 In fact, f3 is proportional to the inverse of the Hawking temperature T. 
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anisotropic order parameter is small and Yang-Mills coupling is very large. However, the ar- 
gument is quite general. Therefore, for any temperature, we expect that the shear viscosity 
behaves as EqflTJ in IR region. 

This note is organized as follows. In SecJ21 we briefly review Einstein- 577(2) Yang-Mills 
system and its analytic solutions. This analytic solutions will be used to provide an example 
for our claim in Sec J3. 11 In Sec|3j we will derive RG equation for the shear viscosities in the 
symmetry broken direction and explore their properties. 



2 Analytic Solutions in Einstein-577(2) Yang-Mills Sys- 
tem and Non-universal Values of the Shear Viscosity 

In this section, we briefly review the analytic solutions in Einstein- <S77(2)-Yang-Mills, the 
dual gravity system of anisotropic superfluids. We mostly follow [19] for the discussion and 
this solution will be used to provide a concrete example of non-trivial RG flows of the shear 
viscosities in this system. 



2.1 Holographic Setup and Analytic Solutions in Large Coupling 
Expansion 

The authors in [TjJ] consider Einstein-S'^7 (2) Yang-Mills system of which space-time is asymp- 
totically AdS$. The action is 



S = j d^xV^G (±(R + H) _ J—F^ N F aMN 



(2) 



where M, N... are 5-dimensional space-time indices, a., are SU(2) indices and g is Yang-Mills 
coupling. For further discussion, we choose L = 1. Yang-Mills field strength F MN is given by 

Fmn = 3mA% — djs[A a M — e abc A b M A c N , (3) 

where e abc is anti-symmetric tensor with e 123 = 1. The equations of motion from the action are 
obtained as 

Wmn = Rmn + ^Gmn — K l \ Tmn ~ -TpGMN^j = 0, (4) 
Y aN = V M F aMN - e abc A b M F cMN = 0, (5) 
where Tmn is the energy- momentum tensor, of which form is 

Tmn = (^F MP F^ a — -F P Q a F PQa G M N S j ■ (6) 

The ansatz for the metric and Yang-Mills field are given by 

A = <p(r)r 3 dt + uj(r)r 1 dx, (7) 

dr^ 

ds 2 = -N{r)a 2 {r)dt 2 + + r 2 f-\r)dx 2 + r 2 f 2 {r) (dy 2 + dz 2 ) , 

N(r) 



where r a = 4^ and s a are Pauli-matrices. We never discuss the detailed forms of equations of 
motion here but just provide their solutions 0. 

A known exact solution of the equations of motion is asymptotically AdS^, charged-black- 
brane solution, of which forms are 



2 

Ml-^f), w(r) = 0, (8) 

2 m 2fl 2 a 2 rj L 



a(r) = f(r) = 1 and N(r) = N (r) = r 2 + 



r 2 3r 4 



.,2 _ 4 „u m - -4 _,_ V" 2 ^ 



where fl is chemical potential, is the black brane horizon, a — and m = r\ + 
In the infinite Yang-Mills coupling limit as g — > oo, the last term in N(r) vanishes and the 
solution becomes uncharged. 

To explore this system near critical point, fl = Jx c = 4^, by analytic method , the authors 
in [19] develop double expansion to the metric fields and Yang-Mills field order by order in eD\ 
and a 2 , e is dimensionless small parameter and Di is the 50(3) rotational symmetry breaking 
order parameter, appearing in the non-trivial zero mode of u(r), 



^) = ^7^y+0(e 2 Dl). (9) 



It is convenient to choose a convention that the horizon of the black brane is located at r = 1 
by scaling that r — > r h r, {t,x,y,z} —> — {t, x, y, z} and defining a new chemical potential 
/i = — . The equations of motion enjoy certain scaling symmetry|16[ H~5] . By these, we can 
choose the asymptotic values of cr(r = oo) = 1 and f(r = oo) = 1 on the large r boundary for 
the space-time to become asymptotically AdS^,. The value of chemical potential is taken to be 
/j, — 4 for the dual boundary field theory system to be at the critical point. 

To obtain the corrections in this double expansion, any appearing fields, a(r) in the 
ansatz([7|) can be expanded as 

a(r) = a (r) + eai(r) + e 2 a 2 (r)... (10) 

Each di in the above expression can also be expanded as 

CLi(r) = a it0 (r) + a 2 a i)2 + a A a iA {r)... (11) 

The zeroth order solutions in e is given in Eq®, where only Nq contains the subleading 
correction of a 2 in the sense of the above expansion. iVo^ = y — p-) and N 0ti = for 
% = 4,6.... It turns out that the non-trivial leading order correction is 0(e 2 a 2 ). The leading 
corrections for the metric are given by 

«) - l -^w£*r' 'w^- v |w (12) 



5 For the precise equations of motion, see Sec. 2 in 
6 For numerical approaches, see [16j [17] 
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2 D\ ( (l + 2r 2 ) 1 281 



4 V3r 2 (l+r 2 ) 3 8 1680 



l-^)+0(£ 3 ). (14) 



(13) 



The black brane temperature is changed by the leading corrections as 




(15) 



where T r = - (l — ¥r 

L IT \ i 



) is the critical temperature. The black brane entropy is given by 




rescaled coordinate. 

2.2 Non-Universality in Anisotropic Background 

In anisotropic background developed in the previous subsection, it is manifest that SO (3) 
isometry in the background metric is broken down to SO (2). As long as we are looking at 
a solution with SO (3) symmetry, the universality of the ratio of entropy density and shear 
viscosity holds. This is because shear viscosities rely on the gravitational perturbations in 
tensor modes of SO (3) in the dual gravity. Each tensor mode satisfies a massless scalar field 
equation decoupled from one another, which can ensure the universality. The universality of 
this ratio is lost in the symmetry broken phase. The reason is that once 50(3) symmetry 
is broken into 50(2), the gravitational wave modes in the symmetry broken direction are no 
longer tensor modes in residual 50(2). They will not be decoupled from other fields and in 
fact interact with Yang-Mills fields. In [19J, the authors have computed a deviation from the 
universal value of the shear viscosity in the symmetry broken direction as 



7 Any subleading corrections of Yang-Mills field in a 2 would not contribute to the leading back reactions 
to the metric. The aim of the calculation in [19] is to get metric corrections up to non-trivial leading order 
corrections, 0(a 2 e 2 ). Therefore, it is enough that the Yang-Mills field solutions is evaluated up to 4>i,o and u>i,o 
only. 




(17) 



using double expansion in eD\ and a 2 together with small frequency expansion %. This result 
is obtained via Kubo formula using the holographic computation of retarded Green's function 
from the perturbations, Hmn an d $A M defined as 

Gmn — G MN + huN and A a M = A^p + 6A M , (18) 

where G M \ and A M are background metric and Yang-Mills fields respectively. For example, 
to compute rj xy , one consider h xy together with 5 Ay and 5 Ay (This is minimal set of fields 
interacting one another classified by residual SO (2) and Z2 symmetry in the symmetry broken 
phase). Their equations of motion are given by 

= g»(r) + (- + — + ^ - mV2(r) 1 (r) + (19) 
vK } \r N(r) f{r) 3N{r)a 2 {r) J vK } N 2 {r)a 2 {r) 1 } 



+ 



2a 2 ( y _ u{r)^{r)5A\{r) ivu(r)<P(r)5A 2 y (r) 

r 2 / 2 (r) V y N 2 (r)a 2 (r) N 2 {r)a 2 {r) 



- f(r)u'(r)%(r) 



r f(r) N(r) a(r) J yKJ \N 2 (r)a 2 (r) 
2iu(j)(r)5A 2 (r) 



N 2 (r)a 2 (r) ' 



l2ff ,, , ,1 2/'(r) , iV'(r) , a'(r)\ x/l2/ , , , / v 2 + 2 (r) 



where ^(r) = ^2^) and v is frequency of the fields. The equations are evaluated in frequency 

space and spatial momenta of the fields are turned off, k = 0. The precise forms of the solutions 
of the weak fields are given in |T9j 0. 



3 i?G-Flows of the Shear Viscosity in Einstein- SU (2) Yang- 
Mills System 

3.1 RG Flow of the Shear Viscosity in Symmetry Broken Direction 

As discussed in [HI EES], retarded Green's function for the gravitational perturbation, ^f u (r) = 
h y x on r = r hyper surface is given by the ratio of its canonical momentum 11^ (r) to itself, 

8 The shear viscosity in the direction of residual symmetry is still universal, rj yz = which depends on 

tensor modes in residual SO(2). 

9 The solutions are quite complicated, so we would not show all these here. For detailed solutions, see Eq(34) 
in Sec3.2 and Appendix. C in [19 . In Sec l3.2l in this note, we need to introduce an 0(1) integral constant, A^\, 
which appears in the solution of 5A y {r). This will be determined by boundary condition in 0{v 2 ), but it is 
difficult to determine that analytically. So we leave this as undetermined in the following discussion. In fact, 
it is not relevant for the shear viscosity computation neither at the horizon or on the AdS boundary. 



6 



which is given by 



n*(f) r 3 a(r)N(r)f e (r)d f V u (r) 



One can define a quantity, 77^ on the r = r hyper surface as 



Vxy{r, u) 



IV 



(22) 



(23) 



for the future convenience. Kubo formula for the shear viscosity cares the imaginary part of 
the retarded Green's function only Therefore, the shear viscosity rj xy is given by 



r] xy = Re[fj. 



xy\- 



(24) 



The holographic RG flow equation of the shear viscosity can be derived from the bulk 
equations of motion. Using Eq (l22|) and Eq(l23l. we switch ^f v (r) in Eq( IT9|) to fj xy , then we get 



df1] X y 



IV 



a(f)N(f) 



+ a 2 A(r) 



(25) 



where 



AM ^ rJ W^ Ur)SA-(r) - "'ffffffi + ) . (26) 



N 2 (r)a 2 (r) 



N 2 (r)a 2 (r) 



On the right hand side of Eq( l25i) . two terms in the square bracket will be the only terms re- 
maining when the the black brane stays above its critical temperature, T > T c (The background 
geometry becomes isotropic). In zero frequency limit, v — > 0, these terms vanish. Then, Eq (l25|) 
becomes dffj xy = and the Z?G-flow is trivial. In this case, the shear viscosity becomes the 
universal value for any r = r hyper surfaces. 

However, below the critical temperature, the background geometry becomes anisotropic 
and the A(r) term in Eq fj25j) appears, which comes from interactions with perturbative Yang- 
Mills fields. Therefore, it is manifest that in the zero frequency limit, the holographic RG flow 
equation of the shear viscosity is not trivial as the usual cases [61 IT3] . in fact it flows from its 
UV boundary value ffTT|) to another one all the way to the IR region. 

The shear viscosities deeply in IR region in the dual field theories correspond to those at 
the black brane horizon in the bulk gravity system. In general ground, it is expected (e.g. see 
|31j ) that the shear viscosities at the horizon in SO(3) rotational symmetry broken direction 
will be modified from its universal value by the ratios of the metric components in the unbroken 
symmetry directions to those of the broken directions. In our case, these might be given by 



Vxy{l 



1 Gyy(r 



2k 2 G xx (r = 1) 



and r] x 



1 G, 



1) 



2n 2 G xx (r = iy 



(27) 



1 is the location of the black brane horizon and ^ is the universal value of 



where again r 

the shear viscosity. In the following, we will argue that this expectation is indeed right even 
in the case that the shear viscosity runs non-trivially. 
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We firstly investigate the near horizon limit of the equations of motion f fT9|) . It turns out 
that the contributions from interaction term A(r) becomes relatively weak near horizon. Once 
we impose ingoing boundary condition at the horizon for the fields appearing in Eqf ll9j) . Eq (|20p 
and Eqf liZTj) as \I/ ~ 5A y ~ 5A y ~ (r — l)~ tul3 , where /3 is a positive real number, the terms 
involving SA y and 8A y in Eq fll9l) become less singular than the terms which are proportional 
to ty v and its derivatives. Then, dominant behavior in Eq( IT9|) presents approximately near 
horizon limit of the equation in isotropic phase as 

^ K{r) + ^ + J^L> . (28 ) 
r — 1 i\ z {rja z (r) 

In fact, the solution of the gravitational field, ^ v {r) has the following near horizon form in 
small frequency limit: 

oo 

tt„(r) = A Q (r, u)(r - ^ z/^„(r), (29) 

n=0 

where iVp = iwv-fi^i ' °"o = °"( r = 1) ; an d A)(?*> z/ ) an d ^n( r ) are some regular functions in 
r and their derivatives are regular too. Once we plug Eqf l29p into the definition of retarded 
Green's function fl2"2"|) . it becomes 



r 3 a(r)N(r)f(r) 
7T 

<5 



Gy(r, V) = —2 d r 



/// 



N <70 



ln(r - 1) + lnA (r, v) + In I ^ ^n{r) 



v ra=0 



(30) 



The last two terms in the square bracket are regular functions in r and N(r) presents single 
zero at the horizon. Therefore, the first term in the square bracket is only surviving at the 
horizon, which gives 

G*(r = l,u) = IK ) and Vxy = ^f(r = 1). (31) 



Gyy(r) 



' G xx (r) 



as 

r=l 



The factor / (r = 1) in the Green's function and r\ xy is precisely the ratio 
expected in ( 1271) . 

One can derive the same result by using Eddington-Finkelstein coordinate in the gravity 
system. Requesting ingoing boundary condition is imposing regularity at the black brane 
horizon. Therefore, the any fields at the horizon depends on t and r only through their non- 
singular combinations. The ingoing null coordinate v is given by 

dv 

dv = dt+ (32) 
N(r)a(r) 

and this implies 

- ii/¥„(r) = N(r)a(r)d r ^(r). (33) 
Plugging this relation into f|22|) . one obtains the same Green's function with fl3TT) . 
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3.2 RG-F\ow Example of the Shear Viscosity from Analytic Solu- 
tions 

In this subsection, we provide an example of RG flow from the analytic solutions to support the 
argument in Sec l3.11 The analytic solution of the bulk equations of motion is already obtained 
in [ToJ^I and we will use its result. Up to leading order corrections in the double expansion of 
eDi and a 2 , we obtained a radial flow of the retarded Green's function by using Eq (l2"2l as a 
function of r: 



Gv(r,v) = -7^2 + 



iv a 2 e 2 



2 K/t^ Ki't^ 



56(1 + r 2 ) 4 + 5376(1 + r 2 ) 4 I ' v ~' § 



10752r 2 (r 2 - l) 2 ^ (34) 



+ iDi(1193 - 2380r 2 + 726r 4 - 348r 6 - 87r 8 + 336r 2 (r 2 - l) 2 /n(l + — )) 

+ higher order in u, e or a 2 , 

where ^ = Zzm r ^ 00 \I / ;/ (r), A{q is an integration constant, which is hard to be determined 
analytically!"! So we keep A^ as unknown in this discussion. We note that e 2 o? = 1260 1 7r 7 ( ?^ Tc \ 
so the temperature dependence of the retarded Green's function is encoded in this term. Only 
imaginary part of the retarded Green's function contributes to the shear viscosity, so the first 
term in the square bracket does not contribute since it is purely real. The shear viscosity on 
r = f hyper surface is given by 



1 a 2 e 2 D 1 

T]xy 



( 



2k 2 5376^ 2 (l + f 2 ) 4 I 



10752f 2 (f 2 - l) 2 im 



(35) 



1 

2 



+ 1)1(1193 - 2380r 2 + 726r 4 - 348r 6 - 87r 8 + 336r 2 (r 2 - l) 2 /n(l + — )) 



+ higher order in u, e or a . 

Since the entropy density is constant in radial direction, s = =£, the ratio of the shear viscosity 
and entropy density runs non-trivially too. As r — > 00, the Green's function and the shear 
viscosity become 

G.(r = oo, „) = (l + ^eWbl) and , m (r = oo) = ^ (l + JLeVflj) , (36) 

up to leading order correction. At the black brane horizon f = 1, the retarded Green's function 
and the shear viscosity flows into the values as 

G»(f = l,u) = -^(l + ^e 2 a 2 Dj) and Vxy (r = 1) = JL + ^s 2 a 2 Df ) . (37) 



10 See Appendix.C in Q]5] 

n See Eq(72), Appendix.C in [19] 



9 



The factor causing deviations from the universal value of the shear viscosity at the horizon is 

1 + yWDl = f(r = 1) = (38) 
in this perturbative regime, which is exactly as expected in ( 1271) . 
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